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SUMMARY 

Five  alloys  of  engineering  interest,  the  cobalt -base  HS-25  and  H-8187  and  the  iron- 
base  SICROMO-9M,  AM-350,  and  AM-355,  were  corrosion  tested  in  boiling  mercury  in 
the  temperature  range  of  1000°  to  1300°  F  (811°  to  977°  K)  for  times  as  long  as  5000 
hours.  Metallographic,  chemical,  and  physical  analyses  were  used  to  determine  the  ex¬ 
tent  and  the  nature  of  the  corrosion. 

Corrosion  mechanism  and  reaction  kinetics  were  inferred  through  a  comparison  of 
experimental  and  theoretical  kinetic  constants.  The  experimental  kinetic  constants  (i.e. , 
the  time-law  constant  br  and  the  activation  energy  AHr)  were  obtained  from  the  test 
data  by  multiple  regression  analysis.  The  theoretical  kinetic  constants  b^  and  AE^ 
for  each  of  i  steps  of  proposed  corrosion  models  were  obtained  from  transport  theory 
relations  and  the  pertinent  boundary  conditions. 

Changes  of  mechanism  during  the  progress  of  corrosion  were  indicated  for  all  the 
alloys  tested  except  the  AM-350  and  AM-355  alloys.  The  shifts  from  one  corrosion  re¬ 
gime  to  another  were  associated  with  changes  in  the  nature  and  length  of  the  corrosion- 
affected  zone.  Three  corrosion  regimes  were  characterized,  linear  I,  parabolic,  and 
linear  n.  The  probable  rate -determining  step  in  the  two  linear  rate  regimes  is 
boundary-layer  diffusion.  The  probable  rate -determining  step  in  the  parabolic  rate  re¬ 
gime  is  liquid  diffusion. 


INTRODUCTION 

|"™The  present  interest  in  space  power  systems  that  employ  mercury  as  the  thermo-  . 
dynamic  working  fluid  at  temperatures  in  the  range  of  1000°  to  1200°  F  (811°  to  922°_K)J 
(refs.  1  and  2)jhas  prompted  the  need  for  an  understanding  of  the  mechanism  and  kinetics 
of  mercury  corrosion.  Such  knowledge  could  be  an  invaluable  guide  in  establishing  oper-  f 


ational  limits,  extrapolating  performance,  and  translating  small-scale  test  results  to 


larger  system  applications  ,£F 

According  to  the  prevailing  concepts  of  reaction  mechanisms,  a  reaction  such  as  the 
corrosion  of  a  solid  material  by  mercury  occurs  as  a  sequence  of  steps.  The  combined 
effect  of  these  steps  determines  the  net  rate  of  the  reaction.  However,  if,  as  is  common, 
one  step  is  significantly  slower  than  the  others,  the  rate  of  the  entire  reaction  is  gov¬ 
erned  by  the  slowest  step.  In  cases  where  the  slow  step  involves  a  diffusion  process 
(e.  g. ,  volume  diffusion  in  a  solid),  the  reaction  is  termed  diffusion  controlled  and  is 
governed  by  the  laws  of  diffusion  kinetics.  Alternately,  if  a  boundary-region  process 
(e.  g. ,  solution)  constitutes  the  slow  step,  the  rate  of  the  overall  reaction  is  determined 
by  the  kinetics  of  this  process.|j 


The  kinetics  of  each  step  of  a  corrosion  mass -transport  reaction  can  be  expressed 
by  the  empirical  relation 


(1) 


where  w.  is  a  measure  of  the  mass  transported,  b.  is  the  time-law  constant  for  the 
th  ^  ^ 

i  step,  and  t  is  time.  The  quantity  k.  is  known  as  the  specific  rate  constant  for  the 
i  step  and  is  a  function  of  temperature  only;  that  is, 


k^  =  ai  exp 


(2) 


where  a.  is  a  constant,  R  is  the  gas  constant,  T  is  absolute  temperature,  and  AH.  is 

1  th  1 

the  "energy  of  activation"  of  the  i  step. 

Therefore,  when  a  corrosion  reaction  is  governed  by  one  controlling  step,  the  mass 
transport  and  the  specific  rate  constant  for  the  entire  reaction  will  be  equal  to,  respec¬ 
tively,  the  mass  transport  w  and  the  specific  rate  constant  k  of  the  rate -controlling 

s  s 

step.  This  relation,  then,  forms  a  basis  for  determining  the  mechanism  and  kinetics  of 
corrosion. 

Only  two  mechanism  studies  of  high-temperature  dynamic  mercury  corrosion  have 
been  reported  (refs.  3  and  4).  In  each  study  the  same  data  were  used,  namely,  the  data 
for  the  corrosion  of  low-carbon  steel  in  mercury  over  the  temperature  range  of  900°  to 
1200°  F  (755°  to  922°  K)  presented  in  reference  5.  Nevertheless,  each  study  reached  a 
different  conclusion  as  to  the  mechanism  of  attack.  The  basic  approach  followed  can  be 
summarized  as  follows: 

(1)  An  assumption  was  made  of  the  rate -controlling  step. 

(2)  A  heat-transfer  -  mass-transfer  analog  equation  was  derived  for  k  . 

s 
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(3)  Estimates  were  made  of  the  mass-transfer  coefficients  and  empirical  constants 
of  the  heat-transfer  -  mass-transfer  analog  equation  and  a  solution  was  obtained  for  kg. 

(4)  The  experimentally  observed  corrosion  rate  was  compared  with  k  for  agree- 
ment. 

(5)  The  activation  energy  derived  from  the  experimental  data  was  compared  for 
agreement  with  the  activation  energy  estimated  to  be  associated  with  the  assumed  rate¬ 
controlling  step. 

In  the  first  study  (ref.  3),  liquid-boundary-layer  diffusion  was  assumed  to  be  the 
rate-controlling  process.  Order -of -magnitude  agreement  was  obtained  between  the  rate 
predicted  by  the  derived  heat -transfer  -  mass-transfer  analog  equation  and  the  observed 
corrosion  rate.  The  order-of -magnitude  agreement  between  observation  and  theory  was 
considered  to  justify  the  general  validity  of  the  premises  made. 

A  second  study  (ref.  4)  used  the  same  data  and  mass-transfer  equation  as  the  first 
study.  It  was  concluded  therein  that  the  difference  between  the  estimated  and  experimen¬ 
tal  activation  energies  and  the  disparity  between  the  observed  and  theoretical  corrosion 
rate  suggested  that  a  process  other  than  diffusion-control  solution  attack  (probably  inter¬ 
granular  attack)  is  the  controlling  factor  in  the  corrosion  of  steel  by  mercury. 

For  completeness,  corroboration  of  a  corrosion  mechanism  by  this  approach  re¬ 
quires  the  formulation  of  a  number  of  corrosion-rate  equations  equal  to  the  total  number 
of  possible  rate -determining  steps.  Concomitantly,  values  of  the  pertinent  material- 
transport  coefficients  and  empirical  constants  are  required  to  solve  each  equation.  In 
many  instances,  these  values  have  not  been  determined  or  are  not  known  accurately. 

In  the  present  study,  a  somewhat  different  approach  was  used.  No  assumption  of  a 
rate -controlling  step  was  made.  Rather,  rate  and  rate -controlling  step  were  inferred 
from  a  comparison  of  experimental  and  theoretical  kinetic  constants.  The  major  prob¬ 
lem  of  this  approach  is  the  paucity  of  available  data  needed  in  the  estimate  of  activation 
energy  for  several  pertinent  processes.  In  brief,  the  approach  was  as  follows: 

(1)  Values  of  the  experimental  kinetic  constants  br  and  AHr  were  obtained  from  a 
multiple  linear  regression  fit  of  the  corrosion  data  to  the  estimating  equation 

AHr 

log  Pr  =  a  +  b  log  t  + -  (3) 

r  r  2.303  RT 

where  pf  is  observed  maximum  corrosion  penetration  or  wall  recession  (assumed  to  be 
proportional  to  wr),  t  is  test  time,  T  is  test  temperature,  and  af  is  a  constant. 
Equation  (3)  was  derived  by  the  elimination  of  k  between  equations  (1)  and  (2)  (where  the 
subscript  r  replaces  i). 

(2)  A  corrosion  model  was  constructed  consistent  with  the  general  nature  of  the  ob¬ 
served  corrosion  attack.  For  each  step  of  the  corrosion  model,  an  equation  was  derived 
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(from  considerations  of  transport  theory  and  the  relevant  boundary  conditions)  in  the 
form  of  equation  (1),  which  yielded  bi  directly. 

(3)  Values  of  AH^  were  obtained  from  the  relation 


AHi  = 


R  d(ln  kt) 


(4) 


which  is  the  differential  form  of  equation  (2). 

(4)  The  theoretical  kinetic  constants  AH^  and  the  b.^  for  each  step  in  the  proposed 

corrosion  mechanism  were  compared  with  AH  and  b  ,  respectively.  Agreement  of 

ij.  r  r 

the  theoretical  kinetic  constants  for  the  i1  step  with  the  respective  experimental  kinetic 

th 

constants  was  taken  to  infer  that  the  i  step  is  rate  determining. 

The  method  of  data  analysis  employed  in  this  approach  provides  a  ready  opportunity 
to  test  for  the  possibility  of  mechanism  changes  during  the  process  of  corrosion.  Such 
changes  may  be  expected,  inasmuch  as  corrosion  induces  physical  changes  in  the  test 
material  which  can  lead  to  changes  in  the  material-transport  path  and,  thereby,  in  the 
rate-controlling  step.  This  process  is  analogous  to  certain  oxidation  reactions  where 
changes  in  oxidation  kinetics  accompany  changes  in  the  thickness  or  adherence  of  oxide 
layers  (ref.  6).  The  possibility  of  mechanism  changes  during  the  course  of  corrosion 
was  explored.  First,  the  data  for  each  test  material  were  grouped  on  the  basis  of  the 
physical  appearance  of  the  corroded  area.  Second,  a  multiple  linear  regression  fit  of  the 
grouped  data  to  equation  (3)  was  made.  Third,  a  statistical  goodness-of-fit  comparison 
was  made  between  the  grouped  and  nongrouped  data. 

For  the  experimental  part  of  this  study,  108  reflux  capsule  tests  were  performed  on 
five  materials  of  engineering  interest,  AM-350,  AM-355,  SICROMO-9M,  HS-25,  and 
H-8187.  The  test  temperature  range  was  1000°  to  1300°  F  (811°  to  977°  K)  and  test 
times  were  as  long  as  5000  hours.  Metallographic  analyses  were  performed  on  all  test 
specimens  to  determine  corrosion  penetration.  In  addition,  chemical  and  physical  analy¬ 
ses  were  performed  on  selected  specimens  to  determine  the  composition  of  the  corrosion 
zones,  the  corrosion  deposits,  and  the  bulk  alloys.  (The  authors  express  their  apprecia¬ 
tion  to  Dr.  Leo  F.  Epstein  of  General  Electric  Vallecitos  Laboratory  and  Dr.  Malcolm 
Parkman  of  Varian  Associates,  Palo  Alto,  California  for  helpful  review  and  for  pointing 
out  an  error  in  activation -energy  estimation  in  the  original  manuscript. ) 
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EXPERIMENT 


Procedure 


The  effect  of  mercury  corrosion  on  HS-25,  H-8187,  SICROMO-9M,  AM-350,  and 
AM-355  was  determined  by  means  of  reflux  capsule  tests  wherein  the  capsule  wall  served 
as  the  test  specimen.  Table  I  lists  the  composition  of  the  test  alloys  and  table  II  lists  the 
forms  in  which  they  were  received,  their  prior  history,  and  their  sources. 


TABLE  I.  -  COMPOSITION  OF  TEST  ALLOYS 


Alloy 

Ni 

Co 

Fe 

Cr 

Mn 

Mo 

w 

Other 

wt.  % 

HS-25 

9.9 

51. 1 

2.5 

19.  8 

1.2 

_ 

15.5 

H-8187 

1.2 

64.6 

1. 1 

20.6 

.6 

0.  1 

11.8 

SIC  ROMO  9M 

90  3 

8  6 

5 

0 

AM-350 

4. 1 

Mj 

75. 1 

16.6 

1. 1 

2.9 

_ 

0.2  (N  +  C) 

AM-355 

4.3 

76.2 

15.5 

1.0 

2.8 

_ 

.2  (N+C) 

TABLE  II.  -  FORM,  HISTORY,  AND  SOURCE  OF  TEST  ALLOYS 


Alloy 

Form 

History 

Source 

HS-25 

1/2-in.  (1.  27-cm)  round 
stock,  centerless  ground 

Solution  treated  at 

2250°  F  (1505°  K) 

Haynes  Stellite  Co. 

H-8187 

2-in.  (5. 08-cm)  hot- 
worked  square 

Solution  treated  at 

2250°  F  (1505°  K) 

Haynes  Stellite  Co. 

SIC  ROMO -9M 

1-in.  (2.  54-cm)  round 
stock 

Annealed  at  1650°  F 
(1171°  K),  cooled 

50°  F  (28°  K) 
per  hour 

Timken  Roller 

Bearing  Co. 

AM-350 

9/16-in.  (1.  43-cm)  rod, 
centerless  ground 

Annealed  at  1425°  F 
(1047°  K),  air  cooled 

A  llegheny  -  Ludlum 
Steel  Corp. 

AM-355 

9/16-in.  (1.  43-cm)  rod, 
centerless  ground 

Annealed  at  1025°  F 
(825°  K),  air  cooled 

All  e  gheny  -  Ludlum 
Steel  Corp. 

Figure  1.  -  Reflux  capsule. 


C -58680 


Capsule  preparation.  -  Capsules  (fig.  1)  were  machined  from  the  as-received  stock 
to  the  nominal  dimensions  of  1/2-inch  (1.  27-cm)  outside  diameter,  1^  inches  (4.  45  cm) 
in  length  and  0.  040-inch  (0. 102-cm)  wall  thickness.  The  inner  surfaces  were  lapped  to 
a  mirror-like  finish.  The  capsules  were  then  degreased  and  chemically  cleaned. 

Next,  the  capsules  were  filled  to  approximately  one-third  their  volume  with  film- 
free,  triple -distilled  mercury  and  sealed  by  electron-beam  welding  in  a  vacuum  of  about 

r  _  o  o 

10  torr  (10  N/m  ).  During  the  welding  process,  which  took  less  than  8  seconds,  the 
capsules  were  held  in  a  water-cooled  jacket  to  prevent  vaporization  of  the  mercury. 

Heat  treatment  of  the  capsule  as  a  result  of  the  weld  was  confined  to  the  vicinity  of  the 
weld.  The  sealed  capsules  were  inspected  radiographically  to  check  the  mercury  level 
and  the  degree  of  wetting  of  the  capsule  wall. 

Capsule  testing.  -  The  capsules  were  tested  in  argon.  Heating  of  the  capsules  was 
accomplished  by  means  of  individual  capsule  furnaces  which  supplied  heat  to  the  lower 
two-thirds  of  the  capsule.  The  furnaces  consisted  of  a  Nichrome  wire  winding  on  a  mag¬ 
nesium  silicate  core  surrounded  on  the  sides  and  bottom  with  0.  375  inch  (0.  952  cm)  of 
magnesium  silicate.  The  tops  of  the  capsules  were  exposed  so  that  heat  was  lost  to  the 
environment  by  natural  convection  and  radiation. 

Test  beds  of  10  furnaces  were  employed  (fig.  2).  Each  furnace  was  individually 
controlled  to  within  ±10°  F  (±5.  5°  K)  of  a  preset  temperature  by  its  own  powerstat. 
Chromel-Alumel  thermocouples  spot  welded  to  the  top  and  bottom  capsule  surface  were 
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used  to  measure  temperature.  The  thermocouple  output  was  continuously  recorded  on  a 
multipoint  strip-chart  recorder. 

Tests  were  conducted  at  1000°,  1100°,  1150°,  1200°,  1250°,  and  1300°  F  (811°, 
866°,  894°,  922°,  950°,  and  977°  K);  test  times  were  as  long  as  5000  hours.  The  tem¬ 
perature  indicated  by  the  bottom  thermocouple  was  used  to  represent  the  test  tempera¬ 
ture. 

Post-test  processing.  -  At  the  conclusion  of  the  tests,  the  capsules  were  allowed  to 
cool  to  room  temperature,  pierced  by  drilling,  and  drained  of  the  mercury.  The  cap¬ 
sules  were  cut  lengthwise  along  the  axis,  and  residual  mercury  was  removed  by  vacuum 
evaporation. 

One  section  of  each  capsule  was  mounted  and  examined  metallographically  in  the 
etched  and  unetched  conditions.  Depths  of  corrosion  penetration  were  measured  with  a 
filar  eyepiece.  Specimens  of  particular  interest  were  analyzed  by  means  of  an  electron- 
beam  probe  microanalyzer  to  determine  the  elemental  composition  of  the  corroded  areas. 
Selected  corrosion  deposits  were  analyzed  chemically  to  determine  their  composition. 


... 

■—"I.  J—  T1 


A  Results^ 


(  Test  data.  A  summary  of  test  conditions,  maximum  corrosion-penetration  depths, 
and  appearance  of  the  corrosion-zone  structure  is  given  in  table  ni.  |  /j  ^ 
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Visual  analysis.  -  When  it  was  removed  from  the  capsules  following  a  test,  the 
mercury  had  a  bright,  silvery  luster.  Shortly  after  exposure  to  air,  the  mercury  began 
to  discolor,  and  a  sludge  appeared  on  the  surface.  The  fluidity  of  the  freshly  drained 
mercury  varied  markedly  and  was  dependent  on  material  and  test  conditions.  In  all 
cases,  a  decrease  in  fluidity  and  an  increase  in  sludge  formation  were  associated  with 
an  increase  in  the  severity  of  corrosion  attack  on  the  capsule  wall  in  the  condensing  re¬ 
gion  and  an  increased  amount  of  crystalline  deposit  at  the  liquid-vapor  interface. 

In  general,  freshly  sectioned  capsules  showed  complete  "tinning"  of  the  wall  by  the 
mercury.  Crystalline  deposits  were  found  along  the  capsule  walls  in  the  liquid  section, 
predominantly  at  the  vapor  -  liquid  pool  interface.  Qualitatively,  there  was  a  close 
relation  between  deposit  buildup  and  corrosion  penetration  of  the  upper  wall.  The  typi¬ 
cal  appearance  of  sectioned  capsules  following  the  removal  of  residual  mercury  is  shown 
in  figure  3. 


Material:  SICR0M0-9M 

Test  temperature:  1100°  F  (866°  K) 
Test  period.-  5000  hours 


AM-350 
1300°  F (977°  K) 
2000  hours 


H-8187 

1300°  F (977°  K) 
5000  hours 


HS-25 

1200°  F (922°  K) 
5000  hours 


Figure  3.  -  Typical  appearance  of  sectioned  capsules  after  removal  of  mercury. 


Metallographic  analysis.  -  The  maximum  corrosion  attack  always  occurred  in  the 
upper  portion  of  the  capsules.  The  variation  of  corrosion  penetration  down  the  capsule 
wall  was  wedgelike;  penetration  decreased  in  a  regular  manner  from  the  top  to  the  liquid 
pool  (fig.  4).  In  general,  examination  of  the  corroded  areas  revealed  no  aberrations 
(i.  e. ,  no  localized  region  of  attack  with  penetration  much  greater  or  much  less  than  the 
surrounding  region).  The  exceptions  to  this  were  the  AM-350  and  AM-355  alloys.  In 
some  instances,  occasional  localized  deep  penetrations  were  observed.  These  anom¬ 
alous  penetrations  were  disregarded  in  estimating  maximum  corrosion  penetration. 
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I  appearance  of  mercury-corroded  areas.  Photomicrographs  of  longitudinal  section  of  capsule  wall.  (Capsule  inner  wall  surface  at 
photomicrograph.) 


The  corroded  areas  presented  one  of  four  typical  appearances:  (1)  uniform  wall  re¬ 
cession,  approximately  a  130°  groove  angle  at  surface  grain  boundaries  (e.g. ,  AM-350, 
fig.  5(a)),  (2)  channel-like  perforation,  approximately  a  50°  groove  angle  at  surface  grain 
boundaries  (e.g.;  SICROMO-9M,  fig.  5(b)),  (3)  a  fine  porous  structure  (e.g.,  HS-25  and 
H-8187,  fig.  5(c)),  and  (4)  a  multifractured  porous  structure  (e.g.,  HS-25,  and  H-8187, 
fig.  3).  The  groove  angles  given  are  approximate,  inasmuch  as  grain-edge  rounding  and 
groove  asymmetry  precluded  accurate  angle  measurements. 

All  materials  except  AM-350  and  AM-355  showed  a  progressive  change  of  the  form 
of  the  corrosion-affected  zone  with  increasing  test  time  and  temperature.  The  HS-25 
and  H-8187  corrosion  zones  progressed  from  a  smooth  wall  to  a  porous  structure  and 
finally  to  a  multifractured  porous  structure.  The  SICROMO-9M  corrosion  zones  pro¬ 
gressed  from  a  smooth  wall  to  grain -boundary  penetrations  which  deepened  and  broad¬ 
ened  with  increasing  time  and  temperature.  For  the  SICROMO-9M  alloy,  in  addition  to 
the  gross  changes  in  the  corrosion  zone,  grain  transformations  were  noted  (fig.  5(b)). 
Although  extensive  wall  recession  was  measured,  the  corrosion-affected  zone  of  the 
AM-350  and  AM-355  alloys  showed  little  change  from  the  initial  smooth  wall  condition 
other  than  a  slight  surface  roughening  and  grain  transformation  in  a  narrow  region, 

o  o' 

1x10  inch  (2.  5x10  cm)  wide,  adjacent  to  the  alloy -condensate  interface. 

Composition  of  corrosion  zones,  deposits,  and  bulk  alloys.  -  The  compositions  of 
the  corrosion  zones,  the  deposits,  and  the  bulk  alloys  were  determined  by  one  or  more 
of  the  following  techniques:  electron-beam  microanalysis,  chemical  analysis,  and  X-ray 
diffraction. 

Typical  corrosion  zone  compositions  were  determined  for  the  alloys  SICROMO-9M, 
AM-350,  H-8187,  and  HS-25  by  electron-beam  microanalysis  scans  across  the  appro¬ 
priate  region  of  selected  samples.  The  SICROMO-9M  analysis  (fig.  6(a))  shows  no  effec¬ 
tive  compositional  change  in  the  corrosion  zone  from  that  of  the  bulk  alloy.  The  analysis 
of  AM-350  (fig.  6(b))  shows  a  depletion  of  manganese  (Mn),  chromium  (Cr),  and  nickel 
(Ni)  in  the  region  adjacent  to  the  alloy-condensate  interface.  Depletion  of  Ni,  Cr,  and 
to  a  lesser  extent  cobalt  (Co),  with  a  corresponding  enrichment  of  tungsten  (W),  is  shown 
by  the  analyses  of  H-8187  and  HS-25  (figs.  6(c)  and  (d)).  The  compositions  of  the  corro¬ 
sion  zones  of  HS-25  and  H-8187,  as  indicated  by  the  electron-beam  microanalyses,  are 
listed  in  table  IV.  Since  the  electron-beam  analyses  of  these  two  alloys  account  for  less 
than  85  percent  of  the  material,  other  means  (discussed  in  the  following  paragraphs) 
were  used  to  establish  the  composition  of  the  corrosion  zones. 

X-ray  diffraction  patterns  were  obtained  from  the  corrosion  zones  of  selected  HS-25 
and  H-8187  capsules.  These  patterns  indicate  the  presence  of  the  intermetallic  com¬ 
pounds  COgW  and  CogW.  Patterns  obtained  from  either  the  as-received  material  or  non- 
corroded  areas  of  the  test  capsules  do  not  indicate  the  presence  of  Co^W  or  C02W. 
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Distance  from  inside  edge,  in. 


2  4  6  8 

Distance  from  inside  edge,  cm 

(dl  HS -25  tested  at  1100°  F  (866°  K)  for  300  hours. 

Figure  6.  -  Electron-beam  microanalysis  of  corrosion  zone  and  matrix. 


TABLE  IV.  -  COMPOSITION  OF  BULK  ALLOY,  CORROSION  ZONE, 


AND  CORROSION  DEPOSITS 


Material 

Analysis 

Co 

Fe 

Ni 

Cr 

Mn 

W 

Mo 

N+C 

wt. 

% 

HS-25 

Bulk 

Chemicala 

51.1 

2.5 

9.  9 

19.8 

1.2 

15.5 

... 

... 

Deposits 

Chemicala 

53.7 

2.7 

14.5 

27.0 

1.4 

1.7 

... 

... 

Corrosion  zone 

Electron  beam*3 

35 

1 

0 

1 

... 

44 

... 

... 

Mass  balance0 

44 

— 

— 

--- 

56 

... 

... 

H-8187 

Bulk 

Chemicala 

64.6 

1. 1 

1.2 

20.6 

0.6 

11.8 

0. 1 

... 

Deposits 

Chemicala 

69.0 

1.3 

1.4 

26.7 

.6 

2.4 

.2 

... 

Corrosion  zone 

Electron  beam*3 

23 

0 

0 

2 

... 

60 

... 

... 

Mass  balance0 

44 

— 

— 

— 

— 

56 

... 

... 

SIC  ROMO -9M 

Bulk 

Chemical3, 

— 

90.3 

— 

8.6 

0.5 

---- 

0.6 

... 

Deposits 

Chemical3 

— 

87.5 

— 

7.6 

--- 

---- 

.3 

... 

AM-350 

Bulk 

Chemical3 

— 

75. 1 

4. 1 

16.6 

1.1 

— 

2.9 

0.2 

Deposits 

Chemical3 

— - 

69.2 

5.4 

22.5 

.8 

— 

.4 

... 

aAnalyses  performed  on  deposits  of  capsules  19,  10,  28,  and  23,  respectively, 
for  HS-25,  H-8187,  SICROMO-9M,  and  AM-350. 

^Analyses  performed  on  corrosion  zone  of  capsules  4,  3,  23,  and  18,  respec¬ 
tively,  for  HS-25,  H-8187,  SICROMO-9M,  and  AM-350. 

Calculated  from  composition  of  bulk  and  deposits  and  assumption  of  conserva¬ 
tion  of  mass. 


TABLE  V.  -  CALCULATED  FLUID  PARAMETERS  FOR  REFLUX 
CAPSULE  TESTS 


Capsule  bottom 
temperature 

Mercury  vaporiza¬ 
tion  rate 

Mercury  film 

thickness 

Average  mercury 
film  velocity 

°F 

°K 

in.  3/hr 

m3/hr 

in. 

m 

in.  /sec 

m/sec 

1000 

811 

0.34 

0.51X10"5 

0.  39X10'3 

0.  99X10"5 

0. 18 

0.46X10'2 

1100 

866 

.45 

.67 

.43 

1. 1 

.22 

.56 

1200 

922 

.62 

.92 

.47 

1.2 

.27 

.69 

1300 

977 

.84 

1.4 

.54 

1.4 

.35 

.89 
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Chemical  analyses  were  made  of  the  bulk  alloys  SICROMO-9M,  AM-350,  AM-355, 
H-8187,  and  HS-25  and  of  deposits  from  selected  test  capsules  of  SICROMO-9M,  AM- 
350,  H-8187,  and  HS-25.  The  results  are  given  in  table  IV.  In  all  cases,  except  for 
SICROMO-9M,  there  was  an  enrichment  of  the  more  soluble  elements,  Ni,  Cr,  or  Mn, 
in  the  deposits  as  compared  with  the  bulk  alloy.  In  the  case  of  SICROMO-9M,  the  de¬ 
posits  were  approximately  the  same  composition  as  the  bulk  alloy. 

A  mass  balance  calculation  was  made  of  the  weight  percentage  of  W  and  Co  remain¬ 
ing  in  the  HS-25  and  H-8187  corrosion  zones  by  using  the  elemental  composition  data  for 
the  bulk  alloy  and  the  corrosion  deposits.  The  following  assumptions  were  made:  (1) 
all  material  removed  from  the  corrosion  zone  became  deposits  and  (2)  the  amount  of 
elements  other  than  W  and  Co  in  the  corrosion  zone  were  negligible.  Listed  in  table  IV 
are  the  mass  balance  estimates  of  the  weight  percentages  of  W  and  Co  in  the  HS-25  and 
H-8187  corrosion  zones.  The  Co-W  ratio,  based  on  the  mass  balance  estimates,  is  2.  5 
for  both  materials.  This  ratio  is  consistent  with  the  X-ray  diffraction  results,  which 
indicated  the  presence  of  both  Co^W  and  Co^W  as  the  major  corrosion  zone  constituents. 

Mercury  vaporization  rate  and  condensate  film  properties.  -  Table  V  presents  the 
calculated  values  of  the  mercury  vaporization  rate  and  condensate  film  thickness  and 
velocity  for  capsule  temperatures  of  1000°,  1100°,  120Q°,  and  1300°  F  (811°,  866°, 
922°,  and  977°  K).  Details  of  the  calculation  are  discussed  in  appendix  A. 

With  the  aid  of  an  X-ray  image  amplifier  the  course  of  the  mercury  condensation 
process  was  followed  in  a  typical  capsule  at  1200°  F  (922°  K)  over  a  period  of  580  hours. 
During  the  first  16  hours  of  operation  dropwise  condensation  was  noted  (fig.  7(a)). 

Drops  formed  and  grew  at  preferential  sites  along  the  upper  portion  of  the  capsule  and 
periodically  ran  down  the  wall.  After  16  hours  no  drops  were  observed;  the  upper  por¬ 
tion  of  the  capsule  appeared  to  be  of  more  or  less  uniform  density  as  viewed  through  the 
X-ray  image  amplifier  (fig.  7(b)).  Mercury  condensation  apparently  occurred  uniformly 
near  the  top  of  capsule  and  the  condensate  traveled  down  the  capsule  wall  in  film  flow. 

An  estimation  of  the  film  thicknes,  of  the  mercury  condensate,  based  on  a  comparative 
optical  density  measurement  with  i  known  thickness  of  lead  used  as  the  standard,  indi- 

_  O  o 

cated  the  film  thickness  to  be  less  ;  xan  1x10  inch  (2.  5xl0”°  cm).  This  thickness  is  to 

-4  -S 

be  compared  with  the  calculated  film  thickness,  5x10  inch  (1.  27x10  cm). 


ANALYSIS 

In  this  investigation  an  attempt  was  made  to  identify,  for  each  material  tested,  the 
probable  rate -controlling  step  at  any  time  during  the  progress  of  corrosion.  The  method 
of  identification  of  rate-controlling  step  was  based  on  a  comparison  of  the  experimental 
kinetic  constants  br  and  AHr,  respectively,  with  the  theoretical  kinetic  constants  b^ 
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and  for  each  of  i  steps  of  a  proposed  corrosion  model.  Agreement  of  the  theoret¬ 
ical  kinetic  constants  of  a  particular  step  with  the  respective  experimental  kinetic  con¬ 
stants  was  taken  to  infer  that  that  step  was  rate  determining.  In  the  following  presenta¬ 
tion,  the  experimental  kinetic  constants  for  each  alloy  were  obtained  by  a  regression 
analysis  data -fitting  technique.  The  theoretical  kinetic  constants  associated  with  pro¬ 
posed  corrosion  models  were  derived  from  transport  theory  relations  for  the  pertinent 
boundary  conditions. 

The  possibility  of  mechanism  changeover  (i.  e. ,  change  of  rate -controlling  step) 
during  the  course  of  the  corrosion  tests  is  also  evaluated.  This  possibility  is  assessed 

(1)  from  the  experimental  results  by  means  of  data  grouping  and  statistical  testing  and 

(2)  from  a  theoretical  consideration  of  the  corrosion  models. 


Experimental  Kinetic  Constants 

Criterion  for  corrosion  attack.  -  The  capsule  method  of  corrosion  testing  precludes 
the  use  of  weight  change  as  the  corrosion  criterion.  Therefore,  in  this  study,  maximum 
corrosion  penetration  was  used  as  the  sole  measure  of  corrosion  attack.  Two  assump¬ 
tions  were  made  in  applying  this  criterion: 

(1)  Penetration  is  proportional  to  the  corrosion  weight  change. 

(2)  The  test  material  is  relatively  homogeneous  from  the  point  of  view  of  corrosion 
attack  (i.  e. ,  there  are  no  local  regions  of  high  corrosion  susceptibility  such  as  micro¬ 
cracks  or  soluble  constituent  concentrations). 

The  first  assumption  is  validated,  at  least  in  the  case  of  attack  leading  to  uniform  wall 
recession,  by  the  principle  of  the  conservation  of  mass;  it  also  finds  support  in  the  gen¬ 
eral  observation  of  increasing  corrosion  deposit  with  increasing  penetration.  The  second 
assumption  is  supported  by  the  fact  that  except  for  AM-350  and  AM-355  alloys,  no  abnor¬ 
mal  localized  regions  of  penetration  were  found  in  the  corroded  areas  of  the  test  cap¬ 
sules.  In  the  case  of  AM-350  and  AM-355,  where  occasional  deep  localized  penetrations 
were  observed,  the  anomalous  penetrations  were  disregarded  in  estimating  maximum 
penetration. 

Data-fitting  techniques.  -  Corrosion  data,  as  a  function  of  time  and  temperature, 
are  customarily  fitted  to  the  empirical  relations: 

b 

pr  =  krt  r  (5) 

B 

log  kp  =  ar  +  —  (6) 
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where  pf  is  corrosion  penetration  (or  weight  loss),  T  is  the  test  temperature  on  the 
absolute  scale,  t  is  time,  and  kf,  ar,  br,  and  Br  are  constants.  The  time  law  (i.  e. , 
the  relation  between  pr  and  t  at  constant  temperature  T)  is  termed  linear  when 
br  =  1  and  parabolic  when  bf  =  1/2.  The  constant  Br  in  the  Arrenhius  relation 
(eq.  (6))  is  proportional  to  the  activation  energy  of  the  overall  corrosion  reaction  (viz, 
AHr  =  4.  58X10"3  Br  kcal/mole  (19.  2  Br  J/mole)). 

When  the  test  data  are  fitted  to  equations  (5)  and  (6)  by  simple  linear  regression 
(least- squares  method),  the  determination  of  br  and  Bf  is  laden  with  certain  difficul¬ 
ties.  First,  the  precision  of  determination  of  br  and  Br  is  generally  poor  because 
the  dependent  variable  pf,  as  a  result  of  uncontrolled  or  unknown  experimental  factors, 
invariably  has  a  large  random  error  (even  though  the  independent  variables  t  and  T 
may  have  relatively  small  errors  of  measurement).  Second,  and  perhaps  as  significant, 
simple  linear  regression  analysis  of  the  experimental  data  yields  values  for  bf  and  Br 
which  are  dependent,  respectively,  on  the  values  chosen  for  t  and  T  (i.  e. ,  dependent 
on  the  arrangement  of  the  experimental  observations,  ref.  ,7). 

These  difficulties  can  be  minimized  or  avoided  to  some  extent  by  the  use  of  multiple 
linear  regression  analysis  (ref.  7).  A  single  regression  equation  containing  all  the  inde¬ 
pendent  variables  is  used  to  express  the  dependent  variable.  By  addition  of  independent 
variables  to  the  regression  equation,  it  is  possible  to  account  for  some  of  the  unexplained 
variation  in  the  dependent  variable  obtained  by  simple  linear  regression  analysis.  This 
resulting  increase  in  the  precision  of  estimate  of  the  dependent  variable  p  is  preferred 
when  analyzing  data  which  exhibit  large  scatter,  such  as  in  the  present  case.  Addition¬ 
ally,  multiple  linear  regression  analysis  affords  a  more  exact  basis  for  estimating  bf 
and  Br  than  simple  linear  regression  and  also  serves  to  summarize  the  results  of  a 
large  number  of  observations  in  a  single  convenient  statement. 

It  must  be  emphasized  that  a  statistical  analysis  of  the  data  does  not  miraculously 
transform  "bad"  data  into  "good”  data.  The  purpose  in  this  analysis  is  to  obtain  the 
best  possible  estimate  of  the  dependent  variable  and  the  kinetic  constants.  The  uncer¬ 
tainty  in  the  estimates  must  always  be  kept  in  mind  and  continually  viewed  with  reference 
to  a  pertinent  error  statistic. 

Analysis  of  the  test  data  by  multiple  linear  regression.  -  The  corrosion  data  were 
treated  by  the  technique  of  multiple  linear  regression.  The  Regression  Analysis  Pro¬ 
gram  with  Iterative  Evaluation  of  Relations  (RAPIER)  (an  unpublished  computer  program 
written  by  Bert  Henry,  Steven  M.  Sidik,  and  R.  Lambird  of  Lewis)  was  used.  The  esti¬ 
mating  equation  used  was  obtained  by  combining  equations  (5)  and  (6)  as  follows: 

Br 

log  pr  =  ar  +  br  log  t  +  —  (7) 
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dL  is  degrees  of  freedom. 
'am-350  and  AM-355  data  combined. 


The  units  arbitrarily  selected  for  the  variables  were  as  follows:  pr,  inchxlO"  ;  t,  hours; 
T,  °K. 

The  results  of  the  regression  analyses  are  presented  in  table  VI.  The  data  were 
grouped  in  the  following  categories:  total,  group  1,  and  group  2.  The  latter  two  catego¬ 
ries  are  discussed  in  the  next  section.  The  23  AM-350  and  the  7  AM-355  data  points 
were  treated  as  one  data  set,  since  compositionally  and  structurally  AM-350  and  AM-355 
are  almost  identical. 

It  is  apparent  from  table  VI  that  for  the  total  data  (i.  e. ,  all  the  data  points  for  a 

o 

given  material)  the  values  of  the  coefficient  of  multiple  determination  It  are  large  and 

p  III 

range  from  0.  76  to  0.  94.  The  large  values  for  indicate  that  the  functional  relations 
selected  for  the  variables  were  reasonably  good  choices,  in  that  (r^  x  loo)  percent  of 
the  variation  in  the  dependent  variable  p  can  be  accounted  for  when  the  data  are  fitted  to 
equation  (7). 

Taken  by  themselves,  the  values  obtained  from  the  total  data  for  Br  ranging  from 

-2610  to  -6380,  which  correspond  to  AH  values  of  -12.0  to  -29.  2  kilocalories  per 
4  5  r 

mole  (-5.  02x10  to  -1.  22x10  J/mole),  appear  to  be  reasonable  for  this  coefficient. 
However,  except  for  AM-350,  the  values  obtained  for  bf  do  not  correspond  to  those 
associated  with  the  expected  processes,  namely,  boundary  or  diffusional  processes 
which  have  linear  or  parabolic  time-law  relations,  respectively.  The  br  values  (0.79, 

0.  74,  and  0.  68  for  HS-25,  SICROMO-9M,  and  H-8187)  lie  between  the  expected  singular 
values  of  1  or  0.  5. 


Mechanism  changeover.  -  It  is  probably  more  than  coincidental  that  the  materials 
for  which  anomalous  values  of  br  are  found  also  exhibit  marked  changes  in  the  structure 
of  the  corrosion  zone  with  increasing  time  and  temperature.  Oxidation  studies  indicate 
that  mechanism  changes  occur  as  a  consequence  of  changes  in  oxide-layer  thickness  or 
oxide -layer  adherence.  Analogously,  it  may  be  expected  that  growth  or  alteration  of  the 
corrosion  zone  may  change  the  nature  of  the  material-transport  path  and,  hence,  result 
in  a  change  of  corrosion  mechanism.  H,  in  fact,  the  rate  process  changes  during  the 
course  of  a  test,  for  example  from  linear  to  parabolic,  then  the  value  of  br  derived 
from  the  total  test  data  would  assume  some  average  value  br,  where  1  >  bp  >  0.  5. 

The  hypothesis  that  changes  in  the  thickness  or  structure  of  the  corrosion  zone  are 
accompanied  by  changes  in  corrosion  mechanism  was  tested  as  follows.  First,  the  cor¬ 
rosion  data  for  each  material  were  grouped  according  to  the  corrosion-zone  criterion  as 
described  later  in  this  section.  Second,  a  multiple  linear  regression  analysis  was  per¬ 
formed  on  each  data  group.  Third,  the  grouped-data  analyses  were  compared  to  the 
analyses  previously  performed  on  the  total  data.  For  each  material,  data  exhibiting  the 
initial  corrosion  structure  were  designated  as  group  1  and  those  exhibiting  the  succeed¬ 
ing  structure  as  group  2.  Data  falling  within  the  corrosion  structure  region  of  transition 
were  included  in  both  group  1  and  group  2.  The  grouping  criterion  is  as  follows: 
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(1)  SIC  ROMO -9M 

-3  -3 

(a)  Group  1,  superficial  penetrations,  <1X10  inch  (2.5x10  cm) 

(b)  Group  2,  deep  penetrations,  >1x10"^  inch  (2.  5X10”^  cm) 

(2)  HS-25  and  H-8187  w  - 

(a)  Group  1,  porous  structure 

(b)  Group  2,  multifractured  porous  structure 

The  results  of  the  grouped-data  analyses  for  HS-25,  SICROMO-9M,  and  H-8187  are  given 
in  table  VI. 

A  goodness-of-fit  criterion  is  needed  to  compare  the  grouped-data  analyses  with  the 

o 

total-data  analyses.  The  coefficient  of  multiple  determination  Rm,  a  measure  of  the 
proportion  of  the  variation  in  the  dependent  variable  explained  by  the  regression  equation, 
was  used  to  validate  the  choice  of  functional  relations  for  the  variables.  This  coefficient 
is  a  relative  measure  of  goodness  of  fit.  However,  once  the  form  of  the  estimating 
equation  is  found  to  be  satisfactory,  an  absolute  measure  of  goodness  of  fit  should  be 
used  for  comparison  of  analyses  employing  the  same  estimating  equation.  Therefore, 
the  residual  error  variance  (ref.  7),  a  measure  of  the  closeness  with  which  the  estimated 
values  of  the  dependent  variable  agree  with  the  experimental  values,  is  employed  herein 
as  the  comparison  criterion  for  the  variously  grouped  data. 

o 

An  inspection  of  the  values  of  residual  error  variance  s.  shows  that  grouping 

X* 

the  data  according  to  corrosion  structure  has  in  all  cases  significantly  improved  the  pre¬ 
cision  in  the  determination  of  log  p  .  Although  grouping  of  the  SICROMO-9M  data  into 


_  O  _ 

groups  of  superficial  penetration  (<1x10  in.  (2.  5X10"°  cm))  and  deep  penetration 
-8  -3 

(>1X10  in.  (2.  5x10  cm))  was  arbitrary,  it  fortuitously  proved  to  be  a  good  choice. 

Trial  analyses  performed  on  the  SICROMO-9M  data  (not  shown  in  table  VI)  yielded  a 

-3 

poorer  fit  of  the  data  when  values  either  slightly  greater  or  slightly  less  than  1x10  inch 
were  used  to  define  the  division  between  superficial  and  deep  penetration. 

Best  estimate  of  regression  coefficients.  -  The  grouped  data  for  HS-25,  H-8187, 
and  SICROMO-9M  give  values  of  bf  that  closely  approximate  the  expected  kinetic  values 
of  1.  0  or  0.  5  and  also  yield  reasonable  values  of  AHr-  In  view  of  this  and  since  data 
grouping  can  be  rationalized  on  both  a  structural  and  a  statistical  basis,  the  best  esti¬ 
mates  of  the  coefficients  of  equation  (7)  were  considered  to  be  those  from  the  grouped- 
data  analyses. 

For  AM-350  no  breakdown  of  the  data  into  groups  was  warranted.  In  this  instance, 
for  the  total  data  the  value  of  br  obtained  was  approximately  1.  Further,  no  marked 
changes  were  observed  in  the  AM-350  corrosion  zone  over  the  time  and  temperature 
range  of  the  tests.  The  coefficients  for  AM-350,  therefore,  were  considered  to  be  those 
for  the  total  data  (table  VI). 

For  later  comparison  of  results,  the  best-estimate  values  of  br  and  AHr  and 
their  standard  deviations  are  listed  in  table  VII. 


\-3 
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TABLE  Vn.  -  EXPERIMENTAL  KINETIC  CONSTANTS 


Material 

Data 

Time-law  constant, 

br 

Activation  energy, 
-AHr 

kcal/mole 

J/mole 

HS-25 

Group  1 

0.  58±0. 09 

7.  8±3.  7 

33±16xl03 

Group  2 

. 93±0. 06 

20.  6±3.0 

86±13 

SIC  ROMO -9M 

Group  1 

. 98±0. 08 

25.  0±2.  5 

105±11 

Group  2 

. 59±0.  03 

9. 9±1.  7 

41±7 

H-8187 

Group  1 

. 54±0. 13 

4.  8±3.  7 

20±16 

Group  2 

1.01±0. 19 

17. 1±4. 0 

72±17 

AM-350 

Total 

1.  07±0.  09 

29.  2  ±4.  2 

123±17 

Theoretical  Kinetic  Constants 

'^Proposed  corrosion  models.  -Jtwo  corrosion  models  are  presented.  Each  proposed 
model  fits  one  of  the  two  general  modes  of  corrosion  attack  observed  in  the  test  capsules, 
namely,  uniform  attack  or  selective  attack.  In  these  models,  movement  of  atoms  from 
the  bulk  solid  material  to  the  bulk  liquid  mercury  occurs  through  a  series  of  consecutive 
and  parallel  mater ial-transport  steps. 

XL 

In  the  following  presentation,  S.  is  used  to  designate  the  1 11  step  and  S.  ,  is  used 

XL  X  XL  lj  K 

to  designate  the  k  parallel  path  for  the  iin  step.  Shown  in  figure  8(a)  is  a  schematic 
representation  of  the  transport  steps  proposed  for  the  selective  attack  model: 

(1)  Diffusion  in  the  solid  Sp  consisting  of  the  parallel  steps,  volume  diffusion  Sj 

and  grain-boundary  diffusion  ^ 

(2)  Solution  S2,  consisting  of  the  parallel  steps,  grain  solution  S2  ^  and  grain¬ 

boundary  solution  S2  2 

(3)  Diffusion  in  a  static  liquid  Sg 

(4)  Convective  diffusion  S^ 

Figure  8(b)  illustrates  the  transport  steps  proposed  for  the  uniform  attack  model  (viz, 

Sp  S2,  and  S^  - ^  33 

Generarconditions  and  assumptions.  -  During  a  test,  mercury  vapor  condenses  con¬ 
tinuously  near  the  top  of  the  capsule  (see  fig.  9).  The  condensate  flows  down  the  capsule 
wall  as  a  thin  film  in  laminar  steady-state  flow  and  enters  the  liquid  pool  at  y  =  -l .  It 
was  assumed  that  condensation  at  the  mercury  film  -  mercury  vapor  interface  x  =  h  is 
negligible.  Since  the  measured  temperature  variation  along  the  upper  capsule  wall  is 
small  and  the  heat  flux  is  low,  the  film  temperature  was  assumed  to  be  uniform. 
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X 


Figure  9.  -  Mercury  flow  in  reflux  capsule. 


In  the  region  of  condensation  the  concentration  of  dissolved  material  was  assumed  to 
be  approximately  0  (i.  e. ,  c  «  0  at  y  =  0).  Lastly,  it  was  assumed  that  steady-state 
conditions  prevail  at  any  point  within  the  mercury  film  (i.  e. ,  3c/ 3t  =  0). 

Time-law  constant.  -  A  mass-transport  -  time  relation  for  each  step  of  the  proposed 
models  was  derived  in  the  form  of  equation  (1)  from  a  consideration  of  (1)  Fick's  first 
and  second  laws  of  diffusion,  (2)  the  Navier-Stokes  equation,  (3)  the  theory  of  rate  pro¬ 
cesses,  and  (4)  the  conservation  of  mass.  The  boundary  conditions  for  each  material- 
transport  step  are  the  conditions  which  prevail  when  the  step  under  consideration  is  as- 
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sumed  to  be  rate  controlling.  Only  the  overall  step  was  considered  in  the  derivation  of 
the  time-law  constant  b^.  Parallel  paths  were  not  treated  separately,  since  the  mass- 
transfer  -  time  relation  for  each  path  takes  the  same  form  (see  appendix  B).  A  summary 
of  the  estimated  b ^  values  for  material-transport  steps  Sp  Sg,  S3,  and  S4  is  given  in 
table  VIE. 


TABLE  Vm.  -  THEORETICAL  KINETIC  CONSTANTS 


Material- 

transport 

step 

Time -law 

constant, 

bi 

Iron -base  alloys 

Cobalt -base  alloys 

Activation  energy,  AH^ 

kcal/mole 

J/mole 

kcal/mole 

J/mole 

S1 

0.5 

30 

126  xlO3 

30 

126  X103 

S2 

1.0 

11.0 

46.0 

16.4 

68.5 

S3 

.5 

7.0 

29.3 

9.7 

37.6 

S4 

1.0 

17.3 

72.3 

22.7 

94.9 

Diffusion  in  the  solid  S^:  E  diffusion  through  the  bulk  solid  is  the  rate -controlling 
step,  and  with  the  assumption  that  the  solid  is  a  semi-infinite  medium,  the  boundary 
conditions  are  as  follows: 

for  t  =  0,  c  =  Ck  at  x<0  and  0  >  y  >  -l 

for  t  >  0,  cw  ~  0  at  x  =  0  and  0  >  y  >  -l 

where  c^  and  cw  are  the  concentrations  in  the  bulk  and  at  the  wall,  respectively,  of 
the  diffusing  species.  The  solution  of  this  problem  is 

W1  =  «iV^/2t1/2  (®) 

where  is  a  constant,  Dg  is  the  diffusion  coefficient  for  the  species  in  the  bulk  solid 
(assumed  to  be  independent  of  c).  For  the  selective  attack  model,  if  it  is  assumed  that 
the  penetration  is  proportional  to  Wj,  equation  (8)  may  be  written  as 

Pi  =  ^D^V/2  (9) 

where  p1  is  the  penetration  and  0^  is  the  appropriate  proportionality  constant.  For 
the  uniform  attack  model,  if  it  is  assumed  that  recession  of  the  wall  is  proportional  to 
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the  amount  of  diffusing  substance,  solution  of  the  problem  of  a  moving  boundary  (refs.  8 
and  9)  leads  to  a  relation  identical  to  equation  (9)  with  the  exception  that  the  proportion¬ 
ality  constant  is  different  and  p^  represents  the  amount  of  wall  recession. 

Solution  Sg:  If  solution  is  the  rate -controlling  step,  the  boundary  conditions  for 
t  >  0  are  c  =  c(y)  at  x  =  0  and  0  >  y  >  -l.  The  solute  concentration  c(y)  is  a  function 
of  y  as  determined  by  the  convective  diffusion  parameters  of  the  flowing  mercury  film. 
The  amount  of  a  given  species  that  passes  into  solution  per  unit  time  per  unit  surface 
area  of  solid  is 


G  =  ksr  [co  -  c(y)]  (10) 

where  kgr  is  a  solution  rate  constant  (a  function  only  of  temperature)  and  cQ  is  the 
equilibrium  solubility  concentration  of  the  solute  in  mercury. 

If  the  rate  of  corrosion  penetration  or  wall  recession  is  assumed  to  be  proportional 
to  the  flux  of  atoms  through  the  surface  into  the  mercury,  then  the  surface  flux  is  given 
by 


(ID 


where  is  a  proportionality  constant  and  Pg  is  the  corrosion  penetration  or  wall  re¬ 
cession.  For  steady-state  conditions, 


Equating  equations  (10)  and  (11)  and  integrating  with  respect  to  t  yield  the  following 
kinetic  relation  for  the  solution- controlled  process: 


P2  =  “[co  “  c(y3t 
a2 


(12) 


Diffusion  in  static  liquid  S^:  If  diffusion  along  quiescent  liquid-filled  channels  or 
networks  is  the  rate-controlling  step,  the  boundary  conditions  for  t  >  0  are: 

c  =  cQ  at  x  =  pg  and  0  >  y  >  -l 

c  =  c'(y)  at  x  =  0  and  0  >  y  >  -l 
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where  p^  is  the  depth  of  greatest  penetration  of  the  channel  and  c’(y)  is  the  solute  con¬ 
centration  at  the  channel  mouth  (a  function  of  y  as  determined  by  the  convective  diffu¬ 
sion  parameters  of  the  mercury  film). 

The  flux  of  solute  from  the  head  to  the  mouth  of  the  channel  is  inversely  proportional 
to  the  length  of  the  diffusion  path  and  directly  proportional  to  the  concentration  difference 
as  follows: 


Dz[co  -  c'<y>] 

?3 


(13) 


where  is  the  diffusion  coefficient  for  the  solute  in  mercury  assumed  to  be  indepen¬ 
dent  of  c. 

If  the  penetration  rate  is  assumed  to  be  proportional  to  the  flux,  then 


j  =  a 


3 


d]?3 

dt 


(14) 


where  is  a  proportionality  constant.  Substituting  the  relation  for  j  of  equation  (13) 
into  equation  (14)  and  integrating  with  respect  to  t  result  in  the  kinetic  relation  for  the 
static  liquid-diffusion-controlled  process 


P3  = 


tl/2 


Convective  diffusion  S^:  If  convective  diffusion  is  the  rate -controlling  step,  the 
boundary  conditions  for  t  >  0  are,  effectively, 


and 


c  =  cQ  at  x  =  0  and  0  >  y  >  -l 


—  =  0  at  x  =  h  and  0  >  y  >  -l 
dx 

The  concentration  distribution  for  a  solute  diffusing  through  a  flowing  thin  film  in 
laminar  flow  and  steady  state  is  given  in  reference  10  as 
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(16) 


3 v(-  -  —  =  D, 

\h  2h  /  ^  3x2 

where  v  is  the  average  film  velocity  (eq.  (A3))  and  h  is  the  film  thickness  (eq.  (A4)). 

Since  the  liquid  film  is  flowing  at  relatively  low  velocities  (see  table  V),  a  fully  developed 

diffusion  regime  will  be  established  across  the  thickness  of  the  film  a  short  distance  from 

the  top  of  the  film.  Calculations  indicate  that  the  thickness  of  the  diffusion  boundary  layer 

-3  -2 

becomes  comparable  to  the  film  thickness  at  y  ~  4x10  inch  (10  cm). 

The  convective  diffusion  mass  transfer  problem  of  S^  is  analogous  to  the  problem 
of  heat  transfer  to  a  fluid  in  laminar  flow  along  a  wall  which  is  maintained  at  uniform 
temperature.  The  solution  to  the  heat-transfer  problem  (ref.  11),  when  applied  to  the 
convective  diffusion  problem,  yields  the  following  expression  for  the  flux  of  mass  at  the 
surface: 


D.c  f 
l  o 

h 


e'n« 


(17) 


-2 

for  values  of  £  greater  than  0. 12  (i.e. ,  y  >  10  cm)  where  n  and  f  are  constants 
and  £  =  D^y/3vh2. 

If  it  is  assumed  that  the  penetration  rate  is  proportional  to  the  flux, 


(18) 


where  a 4  is  a  proportionality  constant  and  p4  is  the  penetration.  Substituting  the  ex¬ 
pression  for  the  flux  (eq.  (17))  into  equation  (18)  and  integrating  with  respect  to  time 
yield  the  kinetic  relation  for  the  convective  diffusion  controlled  process 


P4  = 


D;fC 
L  0 

ha4 


e"n^ 


t 


(19) 


Activation  energy.  -  The  activation  energy  AH^  associated  with  each  material- 
transport  step  is  defined  by  equation  (2).  For  purposes  of  evaluating  AH^,  the  differen¬ 
tial  form  of  equation  (2),  namely, 


R  d(ln  kj) 


T 


=  AHi 


(4) 
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was  used.  Expressions  for  k^,  kg,  kg,  and  k^  (multiplied  by  a  proportionality  constant) 
are  obtained  directly  from  equations  (9),  (12),  (15),  and  (19),  respectively.  The  expres¬ 
sions  for  k^,  evaluated  at  y«0,  were  introduced  into  equation  (4)  and  the  differential 
operation  was  performed.  The  AIT  values  obtained  are  as  follows  (see  appendix  C  for 
details  of  the  calculation): 


-AH1  =  0.5Qj 

(20) 

-ah2  =  +  Qi 

(21) 

-AHg  =  0.  5  (Qg  +  Qv  +  Rt) 

(22) 

-AH4  =  Qg  +  Qv  +  2.  75  RT 

(23) 

where  Q ,  is  the  effective  activation  energy  for  chemical  diffusivity  in  alloy  A  (where 

a  a  A 

A  represents  any  one  of  the  four  alloys  tested),  Q  and  QOT,  are,  respectively,  the 

b  bf 

effective  activation  energy  of  solution  and  solution  rate  in  mercury  of  alloy  A,  and 

is  the  activation  energy  for  viscous  flow  in  mercury.  The  average  temperature  over  the 

experimental  range,  900°  K,  is  used  for  T. 

Of  the  four  activation  energies  only  O  can  be  readily  estimated.  The  value  of  Qy 

was  assumed  to  be  1.  3  kilocalories  per  mole  (5. 4x10^  J/mole)  (ref.  12).  Because  of  the 

complex  nature  of  the  diffusion  and  solution  processes  of  multicomponent  systems,  sev- 

A  A  A 

eral  simplifying  assumptions  must  be  made  to  obtain  estimates  of  ,  Q  ,  and  Q  . 

ci  sr  s 

As  a  result,  too  much  confidence  should  not  be  placed  in  the  value  calculated  for  AH.. 

A  1 

The  evaluation  of  was  based  on  the  following  assumptions: 

(1)  Volume  diffusion  S1  ,  is  the  predominant  mode  of  diffusion  under  the  test  con- 
ditions. 

(2)  The  test  alloys  can  be  treated  as  pseudobinary  systems  (i.  e. ,  Co-Cr  or  Fe-Cr 
systems). 

(3)  The  thermodynamic  activity  coefficients  of  the  major  constituents  of  each  alloy 
are  effectively  constant  with  change  in  composition. 

The  activation  energy  for  volume  self -diffusion  of  the  major  constituent  elements  of  the 
test  alloys  (Fe,  Co,  Cr)  have  comparable  values  (ref.  13)  that  fall  within  the  range  of 
60±10  kilocalories  per  mole  (2.  51±0.  42x10^  J/mole).  Therefore,  the  value  of  the  effec¬ 
tive  activation  energy  of  chemical  diffusivity  for  all  the  test  alloys  was  taken  to  be  ap- 

proximately  60  kilocalories  per  mole  (2.  51x10  J/mole). 

A 

The  evaluation  of  Q  was  based  on  the  following  assumptions: 
s 
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(1)  Only  the  base  element  of  any  given  alloy  (i.  e. ,  Fe  or  Co)  need  be  considered  in 
the  solution  process. 

(2)  The  influence  of  mutual  solubility  effects  on  the  activation  energy  of  solution  is 
negligible. 

The  first  assumption  is  supported  by  the  fact  that  the  alloy  base  element  was  the  major 
component  of  the  corrosion  deposits  (see  table  IV).  The  second  assumption  is  supported 
by  the  results  of  a  study  of  the  solubility  of  the  alloys  HS-25,  H-8187,  and  SICROMO-9M 
in  mercury  (ref.  14). 

The  data  from  reference  14  for  the  cobalt  solubility  of  the  alloys  HS-25  and  H-8187 
as  a  function  of  temperature  are  presented  in  figure  10.  Plotted  on  this  same  figure  are 
the  data  from  references  15  and  16  for  elemental  cobalt  solubility.  The  data  from  refer¬ 
ence  14  for  iron  solubility  of  the  alloy  SICROMO-9M  are  shown  in  figure  11.  Also  plot¬ 
ted  on  figure  11  are  the  data  from  references  15  and  17  for  elemental  iron  solubility. 
Although  AM -350  data  are  given  in  reference  14,  these  data  are  not  presented  herein  be- 


Inverse  temperature,  °K~* 


Figure  10.  -  Solubility  of  cobalt  in  mercury.  For  elemental  cobalt,  log  (ppm  Co)  =  3. 8725 
-  (3585. 3/T);  standard  deviation  of  log  (ppm  Co),  0. 1925;  standard  deviation  of  slope,  560. 8. 
For  cobalt  from  HS-25,  log  (ppm  Co)  »  4. 1730  -  (4093. 4/T);  standard  deviation  of  log  (ppm  Co), 
0.2029;  standard  deviation  of  slope,  1099. 
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Figure  11.  -  Solubility  of  iron  in  mercury.  For  elemental  iron,  log  (ppm  Fe) 

■  2.4396  -  (2409. 1/T);  standard  deviation  of  log  (ppm  Fe),  0. 1911;  standard  devi¬ 
ation  of  slope,  453.9.  For  iron  from  SICROMO-9M,  log  (ppm  Fe)  =  0. 9442 
-(1426. 1/T);  standard  deviation  of  log  (ppm  Fe),  0.2826;  standard  deviation  of 
slope,  752. 


cause  of  their  extremely  large  scatter.  Our  estimates  of  the  regression  line  and  rele¬ 
vant  standard  deviations  for  the  alloy  and  element  data  are  shown  on  each  figure.  Since 
there  are  only  four  data  points  for  the  H-8187  alloy,  no  regression  estimate  was  made. 
The  slopes  of  the  regression  lines  shown  in  figures  10  and  11,  when  multiplied  by  the 
factor  4.  58,  yield  the  appropriate  activation  energy  of  solution  in  units  of  kilocalories 
per  mole. 

The  hypothesis  that  the  slope  of  the  solubility  line  for  the  base  element  of  an  alloy  is 
the  same  as  the  slope  of  the  solubility  line  for  the  pure  element  (i.  e. ,  the  second  as¬ 
sumption)  can  be  tested  by  analysis  of  covariance  (ref.  18).  The  results  of  a  comparison 
of  the  slopes  of  the  regression  lines  are  as  follows.  For  the  HS-25  alloy  slope  compared 
to  the  elemental  cobalt  slope,  the  95-percent  confidence  interval  for  the  difference  be- 
tween  slopes  is  0.  51±2.  53x10  .  For  the  SICROMO-9M  alloy  slope  compared  to  the  ele¬ 
mental  iron  slope,  the  95 -percent  confidence  interval  for  the  difference  between  slopes  is 
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0.  98±1.  86x10  .  Since  the  interval  in  both  cases  includes  zero,  the  hypothesis  cannot  be 
rejected. 

For  the  reasons  just  discussed,  the  activation  energy  of  solution  for  the  cobalt-  and 
iron-base  alloys  was  taken  to  be 

Q?S‘25  «  Q?'8187  «  Q?°  =  16.  4  kcal/mole  (68.  5X103  J/mole) 
s  s  s 

Q8M  «  QtM-350  »  Q^e  =  11.  0  kcal/mole  (46.  OxlO3  J/mole) 

No  data  are  available  for  solution  rate  activation  energy  in  mercury.  This  activa¬ 
tion  energy  is  probably  related  to  the  energy  needed  by  the  dissolving  species  to  diffuse 
through  an  impurity  or  reaction-product  barrier  film  at  the  solid-liquid  interface.  In 
the  tests  reported  herein,  there  was  no  deliberate  addition  of  a  mercury  solution  rate 
inhibitor  (e.  g. ,  zirconium  or  titanium,  which  tend  to  form  nitride  surface  films  (ref.  19)). 
Further,  X-ray  observations  showed  good  wetting  of  surfaces  in  contact  with  mercury, 
which  suggests  the  absence  of  an  oxide  or  oil  impurity  film  at  the  solid-mercury  inter- 

A 

face.  It  was  therefore  assumed  that  Q  is  zero. 

sr 

Listed  in  table  VIII  are  the  AH^  values  for  material-transport  steps  Sj,  S2,  Sg, 
and  S4  calculated  from  equations  (20)  to  (23)  and  the  appropriate  QA  values. 

Variation  of  penetration  as  function  of  y.  -  Where  the  concentration  of  dissolved 
material  is  also  a  function  of  y  (i.  e. ,  rate  control  by  S2,  Sg,  or  S4),  a  variation  of 
corrosion  penetration  down  the  capsule  wall  may  be  expected.  Corrosion  penetration 
should  decrease  in  a  regular  manner  (determined  by  the  solution  rate  constants  and  con¬ 
vective  diffusion  parameters)  from  y  =  0  where  c(y)  =  0  to  y  =  -l  where  c(y)  <!  c  . 

Mechanism  changes.  -  The  proposed  corrosion  model  for  selective  attack  implicitly 
presupposes  a  change  of  mechanism  during  the  progress  of  corrosion.  Phenomenologi¬ 
cally,  this  can  be  viewed  as  a  consequence  of  changes  wrought  in  the  length  and  nature  of 
the  material-transport  path.  Analytically,  this  is  made  evident  when  the  kinetic  relation 
for  the  overall  corrosion  reaction  (i.  e. ,  for  all  steps  of  the  model  viewed  together)  is 
examined.  The  overall  kinetic  relation  for  the  selective  attack  model  is 

at  =  yx  p2  +  y2p  +  y3p2  +  y4P  (24) 


transport  steps,  as  described  in  the  section  Time -law  constant,  and  p  is  the  corrosion 
penetration  predicted  by  the  selective  attack  model.  (See  appendix  B  for  derivative  of 
eq.  (24). )  It  is  seen  from  equation  (24)  that  when  p  is  small  either  S2  or  S4  (depend¬ 
ing  on  the  relative  values  of  y2  and  y4)  must  be  the  rate-determining  step.  Depending 
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on  the  relative  values  of  the  constants,  either  or  Sg  will  become  rate  controlling  as 
corrosion  penetration  increases  with  time.  To  recapitulate,  the  selective  attack  model 
predicts  that  the  initial  rate-controlling  step  follows  a  linear  corrosion  rate  (associated 
with  a  boundary -region  process).  Further,  during  the  progress  of  corrosion,  this  initial 
rate-controlling  step  is  succeeded  eventually  by  a  step  that  follows  a  parabolic  corrosion 
rate  (associated  with  a  diffusion  process). 

No  material-transport  path  change  is  postulated  for  the  uniform  attack  model. 
Therefore,  where  the  uniform  attack  model  applies,  no  mechanism  change  may  be  ex¬ 
pected. 


DISCUSSION  Of/rESULTsJ 
Corrosion  Rates  and  Rate-Controlling  Steps 


Inferences  of  the  rate  and  the  rate -controlling  step  in  the  corrosion  attack  of  mer¬ 
cury  on  the  materials  tested  are  presented  herein.  These  inferences  are  based  primar¬ 
ily  on  a  comparison  of  the  experimental  and  theoretical  kinetic  constants  that  are  derived 
in  the  ANALYSIS  section  and  listed  in  tables  VII  and  VIII,  respectively.  Subsidiary  evi¬ 
dence  as  to  the  nature  of  the  rate -controlling  step  is  found  in  a  comparison  of  the  ob¬ 
served  and  predicted  variation  of  corrosion  penetration  along  the  capsule  wall.  The  uni¬ 
versally  observed  wedge-like  penetration  pattern  in  longitudinal  sections  of  corroded  cap¬ 
sules  corresponds  to  the  predicted  variation  of  corrosion  penetration  when  any  step  but 
S^  is  rate  controlling. 

|  Cobalt -base  alloys,  HS-25  and  H-8187.  -  In  the  initially  observed  corrosion  regime, 


data  group  1,  the  corrosion  rate  was  parabolic,  and  the  probable  rate -controlling  step 
was  Sg.  In  the  succeeding  corrosion  regime,! data  group  2,|the~corrosion  rate  was  lin¬ 
ear, 
and  S, 


In  the  case  of  both  alloys,  the  AH  values  fell  between  the  AH-  values  for  S9 

1  1  u 


Because  of  the  relatively  small  difference  between  AH.  for  Sg  and  S^, 


the 


uncertainty  in  the  estimate  of  AHr,  and  the  assumption  made  in  estimating  the  AH^ 
values,  an  unequivocal  choice  of  rate -controlling  step  cannot  be  made. 

Iron-base  alloy,  SICROMO-9M.  -  In  the  initially  observed  corrosion  regime, \data 
group  lj|ftKe  corrosion  rate  was  linear.  The  probable  rate -controlling  step  was  taken  to 
be  S^.  The  judgment  was  somewhat  arbitrary  because,  although  the  AHr  value  was 
closer  to  AH^  than  to  AHg,  it  was  significantly  larger  than  AH^.|  (See  the  discussion 
for  AM-350.  )r5Tthe  succeeding  corrosion  regimejjdata  group  2,’fffTe  corrosion  rate  was 
parabolic,  and  the  probable  rate -controlling  step  was  Sg. 

Iron -base  alloy,  AM -350.  -  Over  the  entire  test  range,  the  corrosion  rate  was  lin¬ 
ear.  The  probable  rate -controlling  step  was  S^.  Again,  as  in  the  case  of  SICROMO-9M,  |\ 
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this  judgment  is  somewhat  arbitrary.  The  result  for  AM -350  and  SIC  ROMO  -  9M-$. 

(AHr  >  AH^  by  8  to  12  kcal/mole  (3x10^  to  5x10^  J/mole))  is  similar  to  thaj| reported  by 
Horsley  (ref.  4j/forjjiihe  corrosion  of  mild  steel  in  mercury  from  900°  to  1200°  f|(755° 
to  922°  K),Jassuming  boundary -layer  controlj/  From  Horsley’s  calculations  AHr  was 
greater  than  AH^  by  11  kilocalories  per  mole  (4.  6x10^  J/mole).  These  discrepancies 
between  experimental  and  theoretical  activation  energy  may  possibly  reflect  inaccuracy 
in  the  iron  solubility  data  or,  conversely,  an  incorrect  corrosion  model. 


Mechanism  Change  and  Corrosion  Regimes 


^Analysis  of  the  experimental  data  indicates  that  mechanism  changes  occur  during 
the  progress  of  corrosion  for  all  alloys  studied,  with  the  exception  of  AM-350.  Theoret¬ 


ical  consideration  of  the  corrosion  models  support  these  results.^  ~ — ^ 

When  a  change  in  mechanism  takes  place,  the  corrosion  reaction  enters  a  new  cor¬ 
rosion  regime.  Three  corrosion  regimes  can  be  characterized  from  this  present  study: 


linear  I,  parabolic,  and  linear  II.  In  the  following  exposition  of  the  corrosion  regimes, 
figure  12  is  used  as  a  graphic  example.  In  this  figure  are  the  1100°  F  (866°  K)  corrosion 
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rate  curves,  derived  from  equation  (7),  and  the  pertinent  values  found  in  table  VI  for 
HS-25,  H-8187,  SICROMO-9M,  and  AM-350.  The  solid  portion  of  each  curve  covers  the 
test  range.  The  dashed  portions  are  conjectural. 

j~~Linear  I.  -  At  the  start  of  each  test  and  for  some  time  thereafter  the  inner  surface 
of  a  capsule  is  relatively  smooth,  with  no  deep  corrosion  penetrations.  As  such,  the 
corrosion  reaction  is  best  represented  by  the  uniform  attack  model.  Accordingly,  the 
initial  corrosion  rate  is  expected  to  be  linear  for  all  materials. 

A  linear  rate  is  maintained  when  the  capsule  wall  exhibits  a  continuous  uniform  re- 
cessionj.  Such  behavior  is  exemplified  by  the  AM-350  alloy,  as  indicated  in  figure  5. 
fjAttack  on  the  AM-350  alloy  progresses  by  grain -boundary  penetration.  Since,  however, 
the  groove  angle  formed  is  approximately  130°,  the  corrosion  attack  moving  on  a  wide- 
angle  front  results  effectively  in  wall  recession,  j 
I  Parabolic.  -  Continuing  corrosion  attack  on  the  materials  tested,  other  than  AM-350, 
results  in  the  formation  of  a  channel-penetrated  or  porous-network  type  of  corrosion 
zone.  This  developing  corrosion  zone  interposes  itself  as  a  new  link  in  the  mass- 
transport  chain.  Transport  of  matter  across  the  channeled  or  porous  region  occurs  by 
diffusion  along  the  quiescent  liquid  paths  created  by  mercury  contained  within  the  chan¬ 
nels  or  the  network.  When  the  corrosion  channels  or  network  reach  a  sufficient  depth, 
the  rate  of  material  transport  along  the  liquid  diffusional  path  becomes  slower  than  along 
other  paths  in  the  transport  chain.  Thus,  the  rate -controlling  step  changes,  and  the  ini¬ 
tially  observed  linear  corrosion  rate  becomes  parabolic.  The  change  of  rate  from  linear 
to  parabolic  is  continuous  since  the  physical  changes  occurring  in  the  corrosion  zone  are 
continuous  and  progressive.!  For  simplicity  of  presentation,  no  attempt  was  made  to  in¬ 
dicate  a  continuous  change  in  rate  on  the  pertinent  curves  in  figure  12. 

JpThe  formation  of  channels  in  the  SICROMO-9M  alloy  can  be  attributed  to  intergranu¬ 
lar  attack  advancing  with  a  small  groove  angle,  approximately  50°.  No  satisfactory 
explanation  can  be  offered  for  the  pronounced  difference  in  magnitude  of  the  groove 
angles  produced  in  the  SICROMO-9M  and  AM-350  alloys,  which  in  the  case  of  the  latter 
alloy  leads  to  wall  recession. 

The  formation  of  a  porous  COgW  and  COgW  network  in  the  HS-25  and  H-8187  alloys 
can  be  attributed  to  leaching  of  the  more  soluble  elements  from  the  corrosion  zone. 

Since  Co^W  and  CogW  are  not  found  in  the  as-received  material  nor  in  the  material  out¬ 
side  the  corrosion  zone  of  the  tested  capsules,  it  can  be  concluded  that  these  compounds 
are  formed  in  situ  by  reaction  of  the  residual  elements  following  the  leaching  process^ 

Figure  12  indicates  that  at  1100°  F  (866°  K)  SICROMO-9M  exhibits  a  linear  corro-  ~~ 
sion  rate  until  about  1400  hours,  after  which  the  rate  becomes  parabolic.  A  comparable 
linear -to -parabolic  rate  change  for  HS-25  and  H-8187  remains  unconfirmed  because  of 
the  lack  of  sufficient  short-time  data.  The  initial  linear  rate  curve  shown  for  these  al¬ 
loys  in  figure  12  is  conjectural. 
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I  Linear  n. 


In  the  case  of  the  HS-25  and  H-8187  alloys,  a  change  in  density  probably 


occurs  within  the  corrosion  zone  as  a  result  of  leaching.  From  mass  balance  considera¬ 
tions]  (see  table  IV),|itcan  be  estimated  that  the  density  of  the  fully  developed  HS-25  and 
H-8187  corrosion  zones  is,  respectively,  0.  24  and  0.  18  the  density  of  the  original  mate¬ 
rial.  The  density  changes  and  the  accompanying  volume  changes  probably  induce 
stresses  within  the  corrosion-zone  network  and  between  the  network  and  the  base  mate¬ 
rial.  When  the  stresses  exceed  the  critical  shear  stress,  cracks  develop  (see  fig.  3, 
p.  10).  I  The  cracks  allow  the  flowing  mercury  direct  access  to  the  base  material.  As  a 
result,  the  diffusional  path  through  the  porous  corrosion  zone  is  bypassed,  and  the  corro¬ 
sion  rate  reverts  from  parabolic  to  linear.  The  cracks,  once  initiated,  can  be  expected 
to  grow  as  corrosion  continues  and  to  be  accompanied  by  spalling;  thus,  the  linear  corro¬ 
sion  rate  is  maintained.  This  particular  corrosion  behavior  is  closely  analogous  to  cer¬ 
tain  oxidation  reactions  where  parabolic-to-linear  rate  changes  occur  when  developed 
stresses  cause  cracking  and  scaling  of  oxide  layer 5  J 

In  figure  12,1jthe  corrosion  rate  of  HS-25  and  H-8187  at  1100°  F  (866°  K)  changes 
from  parabolic  to  linear  at  1100  and  2400  hours,  respectively.  This  change  is  consistent 
with  the  first  appearance  of  cracks  in  the  corrosion  zone  of  test  capsules.  Cracks  ap¬ 
pear  between  1000  and  2000  hours  for  HS-25  and  between  2000  and  5000  hours  for  H-8187 .J1 


($</vyV 


Liquid  Flow  Velocity 


Liquid  flow  velocity  cannot  be  treated  as  an  independent  variable  in  reflux  capsule 
corrosion  tests.  In  these  tests  the  mercury  film  velocity  is  a  function  of  the  vaporization 
rate,  which  in  turn  is  a  function  of  temperature.  Since,  however,  the  temperature - 
velocity  functional  relation  is  relatively  weak  and  the  temperature  range  of  investigation 
is  relatively  narrow,  the  effect  of  velocity  on  the  test  results  was  assumed  to  be  negli¬ 
gible.  This  was  justified  as  follows:  The  data  from  the  linear  rate  group  for  each  mate¬ 
rial  were  fitted  by  multiple  regression  analysis  to  the  estimating  equation 


Br 

log  Pr  =  ar  +  br  log  t  +  —  -  0-  75  log  T 


(25) 


The  estimating  equation  was  derived  from  equation  (19)  with  the  appropriate  substitution 
for  the  constants  D^,  cQ,  and  h  in  terms  of  T  (see  appendix  C).  Equation  (25)  reflects 
the  effect  of  velocity  on  corrosion  with  the  assumption  of  control.  The  residual 
error  variance  and  the  coefficients  br  and  Br  (eq.  (25))  did  not  differ  significantly 
from  the  values  found  in  table  VI  for  the  pertinent  data  fit  to  equation  (7). 

Although  velocity  effects  were  not  discernible  in  this  investigation,  such  effects 
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should  be  manifest  at  high  liquid  flow  velocity.  Speculations  and  observations  on  the  ef¬ 
fect  of  high  flow  velocity  on  erosion  and  corrosion  are  presented  in  the  following  para¬ 
graphs. 

Erosion.  -  Flow  velocities  within  the  reflux  capsules  were  0.  2  to  0.  4  inch  per  sec¬ 
ond  (0.  5  to  1.  0  cm/sec),  too  low  to  incur  erosional  effects.  However,  at  higher  velocity, 
materials  inclined  to  intergranular  penetration  are  likely  to  suffer  erosion  through  grain 
removal.  Further,  as  a  consequence  of  this  erosion,  the  mode  of  corrosion  attack  will 
differ  from  that  observed  in  the  capsules.  In  the  case  of  SICROMO-9M,  for  example, 
initial  intergranular  attack  would  produce  a  layer  of  peripherally  corroded  grains  at  the 
alloy-liquid  interface.  Through  the  action  of  a  high  stream  velocity  these  grains  would 
be  readily  swept  away.  As  a  consequence  of  continuing  erosional  grain  removal,  channel 
formation  would  be  curtailed.  With  channel  formation  curtailed,  mechanism  change  is 
precluded,  and  corrosion  attack  would  continue  in  the  linear  I  regime. 

This  speculation  is  supported  by  the  results  of  SICROMO-9M  forced-convection- 
boiling  mercury  loop  tests  run  at  liquid  flow  velocities  as  great  as  108  inches  per  second 
(274  cm/sec)  (refs.  20  and  21).  No  channeled-corrosion  zone  was  found,  although  inter¬ 
granular  attack  was  evident.  Also,  the  wall  thickness  loss  was  considerably  greater  than 
that  predicted  by  the  capsule  tests  for  comparable  time  and  temperature  conditions. 

Corrosion.  -  The  effect  of  liquid  flow  velocity  on  corrosion  is  related  to  the  nature 
of  the  rate -controlling  step  of  the  corrosion  reaction.  When  boundary- layer  diffusion 
S4  is  rate  controlling,  corrosion  will  be  proportional  to  vn,  where  n  is  1/3  for  laminar 
flow  and  n  is  1/2  to  1  for  turbulent  flow.  When  solution  S2  is  rate  controlling,  equa¬ 
tion  (12)  indicates  that  penetration  is  proportional  to  cQ  -  c(y).  The  term  c(y)  is  depen¬ 
dent  on  solution  rate  and  stream  velocity  (as  defined  by  the  convective  diffusion  param¬ 
eters  of  the  stream)  and  can  range  in  value  from  0  to  c  .  Consequently,  p2  will  vary 
from  a  maximum  when  c(y)  =  0  to  0  when  c(y)  =  c  .  For  two  cases  of  interest,  p2  is 
independent  of  v:  for  reflux  capsules  at  the  point  of  maximum  penetration  (i.  e. ,  y  =  0, 
c(y)  =  0)  and  for  high  flow,  low  solubility  rate  situations  (i.  e. ,  c(y)  =  constant). 

The  results  of  this  investigation  suggest  that,  for  AM-350  and  SICROMO-9M  alloys, 
S4  control  is  likely  in  the  linear  corrosion  regimes;  no  resolution  could  be  reached  on 
the  question  of  S2  or  S4  control  for  the  HS-25  and  H-8187  alloys.  Evidence  supporting 
S4  control  is  found  in  the  results  of  tests  of  forced-convection  corrosion  loops  fabricated 
of  SICROMO-9M  alloy  (refs.  20  and  21)  and  HS-25  alloy  (unpublished  results  of  tests  by 
A.  Vary  of  Lewis).  In  all  these  tests,  mercury  saturation  temperature  was  about 
1100°  F  (866°  K).  The  velocity  range  covered  was  from  12  to  108  inches  per  second 
(30.  5  to  274  cm/sec).  Although  the  penetration  data  as  a  function  of  velocity  exhibit  a 
large  variability,  in  all  cases  corrosion  penetration  is  greater  than  that  expected  for 
linear  corrosion  rates  with  S2  control. 
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<?Q. 


SUMMARY  OF  RESULTS 


|  A  total  of  108  mercury  reflux  capsule  tests  were  performed  on  the  cobalt -base  al¬ 
loys  HS-25  and  H-8187  and  the  iron -base  alloys  SIC  ROMO -9M,  AM-350,  and  AM-355. 
(The  latter  two  alloys  are  treated  as  one  alloy. )  The  test  temperature  range  was  1000° 
to  1300°  F  (811°  to  977°  K),  and  test  times  were  as  long  as  5000  hours. 

Metallographic  analysis  revealed  that  corrosion  resulted  in  either  uniform  attack, 
characterized  by  wall  recession  (AM-350),  or  selective  attack,  characterized  by  a 
leached-out  porous  corrosion  zone  (HS-25  and  H-8187)  or  a  channel-like  grain-boundary- 
penetrated  corrosion  zone  (SICROMO-9M). 

Values  of  the  experimental  kinetic  constants  (i.  e. ,  the  time-law  constant  br  and  the 
activation  energy  AHr)  were  obtained  from  a  multiple  regression  fit  of  the  corrosion 
data  to  the  estimating  equation 


log  p„  =  a„  +  b  log  t  + 


2.  303  RT 


where 


pr  is  the  observed  corrosion  penetration,  t  is  test  time,  T  is  test  temperature, 


R  is  the  gas  constant,  and  ar  is  a  constant. 

The  theoretical  kinetic  constants  b^  and  AHp  for  each  of  i  steps  of  proposed 
corrosion  models,  were  obtained  from  transport  theory  relations  and  the  pertinent 
boundary  conditions. 

Analysis  of  the  experimental  data  indicated  that  mechanism  changes  occurred  during 
the  progress  of  corrosion  for  all  alloys  studied,  with  the  exception  of  AM-350.  Theo¬ 
retical  consideration  of  the  corrosion  models  supported  these  results. 

Inferences  of  rate  and  rate -controlling  step  were  made  by  comparing  AHi  and  bj 
for  each  step  in  the  proposed  corrosion  model  with  AHr  and  br,  respectively.  Because 
of  the  lack  of  data  needed  to  estimate  the  activation  energy  for  several  pertinent  pro¬ 
cesses,  there  is  some  uncertainty  in  the  estimated  AH-  values.  Consequently,  the  in¬ 
ference  of  rate -controlling  step  must  be  considered  as  tentative. 

Two  corrosion  regimes  were  inferred  from  the  HS-25  and  H-8187  alloy  data.  The 
first  was  a  parabolic  rate  regime  with  liquid  diffusion  as  the  probable  rate-controlling 
step.  This  regime  was  succeeded  by  a  linear  rate  regime  with  boundary-layer  diffusion 
as  the  probable  rate -controlling  step.  Although  not  experimentally  verified,  there  is 
good  reason  to  believe  that  the  parabolic  regime  is  preceded  by  a  linear  rate  regime. 

Two  corrosion  regimes  were  inferred  from  the  SICROMO-9M  alloy  data.  The  first 
was  a  linear  rate  regime  with  boundary -layer  diffusion  as  the  probable  rate-controlling  I 
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step.  This  regime  was  succeeded  by  a  parabolic  rate  regime  with  liquid  diffusion  as  the 
probable  rate-determining  step. 

A  single  corrosion  regime  of  linear  rate  with  boundary -layer  diffusion  as  the  prob¬ 
able  rate -controlling  step  was  inferred  in  the  case  of  the  AM-350  alloy. \)  ^ 
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APPENDIX  A 


MERCURY  VAPORIZATION  RATE  AND  CONDENSATE  FILM  PROPERTIES 

A  determination  was  made  of  the  temperature  profile  along  the  capsule  and  the  heat 
loss  from  the  top  of  the  capsule  at  several  temperature  levels  in  order  to  estimate  the 
mercury  vaporization  rate  in  the  reflux  capsule  tests.  A  typical  capsule  was  used  with 
five  thermocouples  spot  welded  to  the  wall  in  the  positions  shown  in  figure  13.  Other 


0  1  2  3  4  5  6 

Distance  from  bottom  of  furnace,  cm 

Figure  13.  -  Typical  temperature  profiles  for  a  capsule  containing  mercury. 
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Capsule  temperature 


thermocouples  were  distributed  throughout  the  capsule  furnace  assembly.  The  setup  was 
operated  at  1000°,  1100°,  1200°,  and  1300°  F  (811°,  866°,  922°,  and  977°  K)  with  the 
capsule  one-third  filled  with  mercury.  The  capsule  was  then  removed,  emptied  of  mer¬ 
cury,  and  resealed  in  vacuum  with  an  electron-beam  welder.  The  capsule  was  re¬ 
installed  in  the  furnace  and  operated  over  a  range  of  temperature  from  900°  to  1600°  F 
(755°  to  1144°  K). 

The  temperature  profiles  for  the  runs  made  with  mercury  in  the  capsule  are  shown 
in  figure  13.  The  slight  temperature  rise  indicated  along  the  lower  two -thirds  of  the 
capsule  is  probably  a  result  of  the  proximity  of  thermocouples  2  and  3  to  the  heater.  The 
temperature  drop  observed  over  the  top  one-third  of  the  capsule  is  most  likely  somewhat 
larger  than  the  true  drop,  because  of  heat -conduction  loss  from  convective  cooling  of  the 
thermocouple  wire. 

Temperature  profiles  obtained  with  the  evacuated  empty  capsule  are  shown  in  fig¬ 
ure  14.  The  profiles  show  a  nearly  uniform  temperature  for  the  lower  portion  of  the 
capsule  and  a  temperature  drop  of  200°  to  400°  F  (1110  to  222°  K)  for  the  upper  portion 
of  the  capsule.  The  anomalous  temperature  drop  indicated  by  the  top  thermocouple  (5) 
probably  results  from  two  factors.  The  first  factor  is  an  added  resistance  to  heat  con¬ 
duction  at  the  thermocouple  location  produced  by  the  dimple  machined  in  the  surface  of 
the  capsule.  The  wall  thickness  at  the  dimple  was  0.  010  to  0.  020  inch  (0.  025  to 
0.  051  cm)  as  compared  with  0.  040  inch  (0. 10  cm)  for  the  nominal  wall  thickness.  The 
second  factor  is  an  enhanced  heat -conduction  loss  down  the  thermocouple  wires  caused 
by  convective  cooling.  It  is  assumed  that  the  true  top  temperature  is  that  value  obtained 
by  straight-line  extrapolation  of  the  temperatures  given  by  thermocouples  3  and  4. 

The  heat  loss  from  the  top  of  the  capsule  was  calculated  from  the  following  equation: 

q  =  kA  (Al) 

L 

where  k  is  the  thermal  conductivity,  A  is  the  cross-sectional  area  of  the  capsule  wall, 
and  AT/L  is  the  slope  of  the  linear  portion  of  the  temperature -profile  curve  between 
the  thermocouples  3  and  5  of  figure  14.  Heat  loss  from  the  capsule  top  is  plotted  in  fig¬ 
ure  15  as  a  function  of  true  top-surface  temperature  of  the  empty  capsule. 

Since  the  heat  loss  from  the  top  of  the  capsule  for  any  given  top  temperature  is  inde¬ 
pendent  of  whether  the  capsule  is  filled  or  empty,  values  of  heat  loss  for  the  filled  cap¬ 
sule  tests  can  be  obtained  directly  from  figure  15. 
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700  800  900 

Temperature  at  top  of  capsule,  °K 

Figure  15.  -  Heat  loss  from  capsule  top  as  function  of  temperature. 


The  vaporization  rate  m  was  obtained  from  the  relation 

m=-9_ 

V 


(A  2) 


where  H  is  the  latent  heat  of  vaporization  and  p  is  the  density  of  mercury.  Calculated 
V  0  0 

values  of  m  are  given  in  table  V  for  capsule  bottom  temperatures  of  1000  ,  1100  , 

1200°,  and  1300°  F  (811°,  866°,  922°,  and  977°  K). 

Average  velocity  v  over  the  film  cross  section  is  (ref.  11,  p.  671) 

v  =  ^  (A3) 

Zv 


where  h  is  the  film  thickness,  v  is  the  kinematic  viscosity,  and  g  is  the  gravitational 
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constant.  A  relation  for  the  film  thickness  can  be  obtained  as  follows.  The  average 
velocity  is  equal  to  the  vaporization  rate  m  divided  by  the  flow  cross-sectional  area  of 
the  film.  Since  h  «  d,  where  d  is  the  capsule  inside  diameter,  the  flow  cross- 
sectional  area  is  approximately  equal  to  7rdh.  Substituting  m/7rdh  for  v  inequation 
(A3)  and  solving  for  h  yield 


(A4) 


The  calculated  film  thickness  and  average  velocity  are  given  in  table  V  for  capsule 
temperatures  of  1000°,  1100°,  1200°,  and  1300°  F  (811°,  866°,  922°,  and  977°  K). 
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APPENDIX  B 


DERIVATION  OF  KINETIC  RELATION  FOR  SELECTIVE  ATTACK  MODEL 

The  kinetic  relation  for  the  selective  attack  model  can  be  obtained  by  considering  the 
overall  corrosion  reaction  in  terms  of  a  flux  of  matter  j  which  has  to  overcome  a  series 
of  resistances  r^.  Each  resistance  is  associated  with  a  corresponding  step  S^.  The 
sum  total  of  all  the  steps  compose  the  mass-transfer  path.  The  driving  force  for  each 
step  is  a  concentration  difference  Ac^.  Thus,  analogous  to  current  flow  in  an  electric 
series  circuit, 


j  = 


(Bl) 


When  parallel  paths  exist,  as  for  example  in  the  case  of  volume  and  grain -boundary 
diffusional  paths  ^  and  g  in  the  solid  diffusion  step  (fig.  8),  the  total  resis¬ 
tance  r  j  is  given  by  the  expression 


rl  rl,l  +  rl,2 

Consequently,  when  r^  g,  ^  «  rj  g/2  or  when  rj  j  >  ri  2’  rl  ~  rl  2'  Thus> 

the  total  resistance  may  be  expressed  in  terms  of  the  smaller  of  the  parallel  path  resis¬ 
tances. 

Two  types  of  resistances  may  be  considered:  (1)  resistance  r^  ^  associated  with 
steps  involving  a  boundary  region  (e.  g. ,  Sg  and  S^)  and  (2)  resistance  r^  ^  associated 
with  steps  involving  diffusional  processes  (e.g. ,  and  Sg).  The  relation  between  r^  ^ 
and  the  depth  of  corrosion  penetration  p  is 

ri,  d  =  yip  <B2> 

i.1. 

where  is  a  characteristic  constant  for  the  i  n  step.  Resistance  r^  ^  may  be  ex¬ 
pressed  simply  as 


(B3) 


When  the  flux  of  matter  is  constant, 
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where  <p  is  a  proportionality  constant.  Combining  equations  (Bl)  and  (B4)  and  making 
the  appropriate  substitution  for  the  ri  terms  yield  the  following  corrosion-rate  expres¬ 
sion  for  the  selective  attack  model: 


£Aci 

,,*  = _ tl _ 

dt  yxp  +  r2  +  r3P  +  y  ^ 


(B5) 


If  it  is  assumed  that  Ac^  is  constant,  integrating  equation  (B5)  and  combining  constants 
into  a  new  constant  a  (the  integration  constant  being  0  since  p  =  0  at  t  =  0)  yield  the 
following  kinetic  relation: 


APPENDIX  C 


ACTIVATION-ENERGY  CALCULATIONS 

The  following  equation  was  used  to  calculate  the  activation  energies  for  each 
material-transport  step: 

R  d(ln  k.) 

AH.  - - L  (4) 

dl 

T 

Listed  below  are  the  pertinent  expressions  for  obtained  from  equations  (9),  (12),  (15), 
and  (19)  and  evaluated  at  the  point  of  maximum  attack.  An  appropriate  proportionality 
constant  e^  is  included: 


ki  = 

1  Is 

(Cl) 

k2  -  e2ksrco 

(C2) 

k3  =  e3®(c0I1/2 

(C3) 

o 

Q 

VLl 

II 

(C4) 

Diffusion  in  Solid  S1 


T 

The  temperature  dependence  of  D  is  as  follows: 

s 

Ds  =  Do  exP 


(C5) 


(C6) 
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where  Dq  is  a  constant  and  is  the  effective  activation  energy  for  diffusivity  in  al¬ 
loy  A.  Substituting  for  Dg  in  equation  (C5)  and  performing  the  differential  operation 
yield 

-AH1  =  0.  5  qJ  (C 


Solution  S< 


AH2  = 


R  d(ln  £2ksrc0) 


The  temperature  dependence  of  k  and  c  is  given  by  the  following  relations: 

“  bl  O 


ksr  =  X  exp 


c  =  A*  exp 


(CIO) 


A  A 

where  A  and  A  are  constants  and  Q  and  Q  are,  respectively,  the  effective  acti- 

&  br 

vation  energy  of  the  solution  and  the  solution  rate  in  mercury  of  the  elements  of  alloy  A. 

Substituting  for  k  and  c  in  equation  (C8)  and  performing  the  differential  operation 
sr  u 


-AH2  =  qA  +  qA 


(Cll) 


Diffusion  in  Static  Liquid  S3 


ah3  = 


_Rdtah(DzC0)1/2] 


(C12) 


The  temperature  dependence  of  can  be  approximated  by  the  Stokes -Einstein  relation: 
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D,  =  erp£?s\  (C13) 

6’Vo  \RT/ 

where  1^  is  the  Boltzmann  constant,  rjQ  is  a  constant,  rQ  is  the  radius  of  the  diffusing 
substance  (approximately  the  same  for  all  metal  atoms  considered  in  this  study),  and  Qy 
is  the  activation  energy  for  viscous  flow.  Combining  equations  (C12),  (CIO),  and  (C13) 
and  differentiating  yield 

-AH3  =  0.  5(Q^  +  +  RT)  (C14) 

Convective  Diffusion  S4 

Rd(ta£4“ir) 

AH,= - X - ?L_Z  (C15) 

4  di 

T 

An  analysis  of  the  mercury  vaporization  data  (table  V)  indicates  that  h  is  proportional 
to  T1’75.  Substituting  the  appropriate  temperature-dependent  relations  for  D^,  cq, 
and  h  into  equation  (C15)  and  differentiating  give 

-AH4  =  Q^  +  Qv+  2.75  RT  (C16) 
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APPENDIX  D 


SYMBOLS 


A  cross-sectional  area  of  capsule 

wall, 

a  constant 

Br  regression  coefficient 

b  time-law  constant 

9 

c  concentration,  kg/m 

c  equilibrium  solubility  concentra- 

0  3 

tion,  kg/m 

D,  diffusion  coefficient  in  liquid, 

i  O 

m^/hr 

D„  diffusion  coefficient  in  solid, 

S  P 

nr /hr 

d  capsule  inside  diameter,  m 

f  constant 

o 

G  solution  mass  flux,  kg/(m  )(hr) 

2  2 

g  gravitational  constant,  (N-m  )/kg 

Hy  latent  heat  of  vaporization,  J/kg 

AH  activation  energy,  J/mole 

h  condensate  film  thickness,  m 

2 

j  mass  flux,  kg/(m  )(hr) 

k  specific  rate  constant 

k^  Boltzmann  constant,  J/°K 

kOT.  solution  rate  constant,  m/hr 

sr 

L  linear  portion  of  temperature  pro¬ 

file,  m 

l  vertical  length  of  condensate 

film,  m 

3 

m  vaporization  rate,  m  /hr 


n  constant 

p  corrosion  penetration  or  wall  re¬ 

cession,  m 
A 

effective  activation  energy  for 
chemical  diffusivity  in  alloy  A, 
J/mole 

A 

Qg  effective  activation  energy  for 

solution  of  elements  of  alloy  A, 
J/mole 

A 

Q  effective  activation  energy  for 

sr 

solution  rate  of  elements  of 
alloy  A,  J/mole 

Qy  activation  energy  for  viscous  flow, 

J/mole 

q  heat  loss,  J/hr 

R  gas  constant,  J/(°K)(mole) 

o 

Rm  coefficient  of  multiple  determina¬ 

tion 

r  mass  flux  resistance,  hr/m 

rQ  radius  of  diffusing  species 

S  mate  rial -transport  step 

s  standard  deviation 

2 

s  variance 

T  absolute  temperature,  °K 

t  time,  hr 

v  liquid  flow  velocity,  m/hr 

v  average  film  flow  velocity,  m/hr 

w  mass  transport,  kg 
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a,P,Y, 1 

constants 

d 

diffusional  process 

e’VxJ 

i 

transport  step 

K 

thermal  conductivity, 

k 

transport  step  parallel  path 

J/(m)(hr)(°K) 

r 

derived  from  corrosion  data 

V 

O 

kinematic  viscosity,  m  /hr 

s 

rate-controlling  transport 

% 

convective  diffusion  parameter 

step 

P 

3 

density,  kg/m 

w 

wall 

Subscripts: 

b 

bulk;  boundary -region  process 

E-2109 
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